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ABSTRACT 

We present the Sunyaev-Zeldovich (SZ) signal-to-richness scaling relation (7500 - A'200) for the MaxBCG cluster catalogue. Employing a multi- 
frequency matched filter on the Planck sky maps, we measure the SZ signal for each cluster by adapting the filter according to weak-lensing 
calibrated mass-richness relations (A'200 - ^soo)- We bin our individual measurements and detect the SZ signal down to the lowest richness systems 
(A'200 = 10) with high significance, achieving a detection of the SZ signal in systems with mass as low as M500 ~ 5 x 10" Mq. 
The observed ^500 - Moo relation is well modeled by a power law over the full richness range. It has a lower normalisation at given A'200 than 
predicted based on X-ray models and published mass-richness relations. An X-ray subsample, however, does conform to the predicted scaling, and 
model predictions do reproduce the relation between our measured bin-average SZ signal and measured bin-average X-ray luminosities. At fixed 
richness, we find an intrinsic dispersion in the F500 - A'200 relation of 60% rising to of order 100% at low richness. 

Thanks to its all-sky coverage, Planck provides observations for more than 13,000 MaxBCG clusters and an unprecedented SZ/optical data set, 
extending the list of known cluster scaling laws to include SZ-optical properties. The data set offers essential clues for models of galaxy formation. 
Moreover, the lower normalisation of the SZ-mass relation implied by the observed SZ-richness scaling has important consequences for cluster 
physics and cosmological studies with SZ clusters. 
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1. Introduction 

Galaxy cluster properties follow simple scaling laws (see e.g. 
Rosati et al. 2002; Voit 2005, for recent reviews). This attests 
to a remarkable consistency in the cluster population and moti- 
vates the use of clusters as cosmological probes. These scaling 
laws also provide important clues to cluster formation, and rela- 
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tions involving optical properties, in particular, help uncover the 
processes driving galaxy evolution. 

The Sunyaev-Zeldovich (SZ) effect (Sunyaev & Zeldovich 
1972; Birkinshaw 1999) opens a fresh perspective on cluster 
scaling laws, and the advent of large-area SZ surveys furnishes 
us with a powerful new tool (Carlstrom et al. 2002). Proportional 
to ICM mass and temperature, the thermal SZ effect probes the 
gas in a manner complementary to X-ray measurements, giving 
a more direct view of the gas mass and energy content. Ground- 
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Fig. 1. Individual scaled SZ signal measurements, Fsoo, for the 
MaxBCG catalogue as a function of richness A^20o- We do not 
plot individual error bars to avoid saturating the figure. The er- 
ror bar drawn in the upper left represents the median uncer- 
tainty over the entire population; in general, the uncertainty in- 
creases towards low richness. The SZ signal measurements are 
expressed as the Compton y parameter integrated over a sphere 
out to 7?5o(), scaled in redshift according to the self-similar model 
and placed at a fiducial angular distance of 500 Mpc. Each point 
represents the result of the matched filter applied to an individ- 
ual cluster in the catalogue. Upward pointing arrows indicate 
values beyond the plotted range. The radius Rsm, and hence the 
filter size, is set from the mass of each cluster determined via the 
weak-lensing calibrated M500 - A^20{) relation given by Johnston 
et al. (2007). The results are nearly the same for the relation 
given by Rozo et al. (2009) 



based instruments, such as the Atacama Cosmology Telescope 
[(ACT), (Swetz et al. 2008)], the South Pole Telescope [SPT, 
(Carlstrom et al. 2009)] and APEX-SZ (Dobbs et al. 2006), are 
harvesting a substantial crop of scientific results. They are pro- 
ducing, for the first time, SZ-selected catalogues and using them 
to constrain cosmological parameters (Staniszewski et al. 2009; 
Marriage et al. 2010; Sehgal et al. 2010; VanderUnde et al. 2010). 

The Planck^ consortium has published its first scientific 
results (Planck Collaboration 2011a) and released the Planck 
Early Release Compact Source Catalogue (ERCSC) (Planck 
Collaboration 2011c), which includes the Planck Early SZ 
(ESZ) all-sky cluster list (Planck Collaboration 201 Id). Planck 
(Tauber et al. (2010); Planck Collaboration (201 la)) is the third 
generation space mission to measure the anisotropy of the cos- 
mic microwave background (CMB). It observes the sky in nine 
frequency bands covering 30-857 GHz with high sensitivity and 
angular resolution from 31 '-5'. The Low Frequency Instrument 
(LFI; Mandolesi et al. (2010); Bersanelli et al. (2010); Mennella 
et al. (2011)) covers the 30, 44, and 70 GHz bands with am- 
plifiers cooled to 20 K. The High Frequency Instrument (HFI; 
Lamarre et al. (2010); Planck HFI Core Team (2011a)) covers 
the 100, 143, 217, 353, 545, and 857 GHz bands with bolome- 
ters cooled to 0.1 K. Polarization is measured in all but the high- 
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est two bands (Leahy et al. (2010); Rosset et al. (2010)). A 
combination of radiative cooling and three mechanical coolers 
produces the temperatures needed for the detectors and optics 
(Planck Collaboration (2011b)). Two Data Processing Centres 
(DPCs) check and calibrate the data and make maps of the sky 
(Planck HFI Core Team (201 lb); Zacchei et al. (201 1)). Planck's 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for galactic and extragalactic astrophysics, 
as well as cosmology. Early astrophysics results are given in 
Planck Collaboration (201 ld)-Planck Collaboration (201 Ih). 

Planck early results on clusters of galaxies are presented 
in this paper and in (Planck Collaboration 201 ld,e,f,g). In the 
present work, we use Planck SZ measurements at the locations 
of MaxBCG clusters (Koester et al. 2007a) to extract the SZ 
signal-richness scaling relation. There are several optical clus- 
ter catalogs (Wen et al. 2009; Hao et al. 2010; Szabo et al. 2010) 
available from the Sloan Digital Sky Survey (York et al. 2000, 
SDSS). For this initial study, we chose the MaxBCG catalogue 
for its large sample size, wide mass range and well-characterized 
selection function, and because its properties have been exten- 
sively studied. In particular, we benefit from weak-lensing mass 
measurements and mass-richness relations (Johnston et al. 2007; 
Mandelbaum et al. 2008a; Sheldon et al. 2009; Rozo et al. 2009). 
A combined SZ-optical study over such a large catalogue is un- 
precedented and Planck is a unique SZ instrument for this task, 
as its all-sky coverage encompasses the complete SDSS area and 
the full MaxBCG cluster sample. 

Our analysis methodology follows that of the accompanying 
paper on the SZ properties of X-ray selected clusters (Planck 
Collaboration 201 If). Although the individual SZ measurements 
in both cases generally have low signal-to-noise, we extract the 
statistical properties of the ICM — mean relations and their dis- 
persion — by averaging over the large sample. The approach 
enables us to study the properties of a much larger and represen- 
tative sample of clusters than otherwise possible. 

The SZ-richness relation adds a new entry to the complement 
of cluster scaling laws and additional constraints on cluster and 
galaxy evolution models. With a mass-richness relation, we can 
also derive the SZ signal-mass relation. This is a central element 
in predictions for the diffuse SZ power spectrum and SZ cluster 
counts. Poor knowledge of the relation represents an important 
source of modeling uncertainty. Low mass systems, for example, 
contribute a large fraction of the SZ power, but we know very 
little about their SZ signal. 

We organise the paper as follows: the next section presents 
the data used, both the Planck maps and the MaxBCG catalogue 
and pertinent characteristics. Section 3 details our SZ measure- 
ments based on a multi-frequency matched filter, and outlines 
some of the systematic checks. In Sect. 4 we present our ba- 
sic results and in Sect. 5 compare them to model expectations. 
Section 6 concludes. 

1.1. Conventions and Notation 

In the following, we adopt a flat fiducial cosmology with 
Q.M = 0.3 with the remainder of the critical density made up by 
a cosmological constant. We express the Hubble parameter at 
redshift z as H{z) = HuE{z) ^{hx 100 kms"' Mpc-')£'(z) with 
h - 0.7. Cluster radii are expressed in terms of /?a, the radius 
inside of which the mean mass overdensity equals A x pdz), 
where pdz) - 3H^{z)/^nG is the critical density at redshift z. 
Similarly, we quote masses as Ma = A(47r/3)7?^pc- We note that, 
in contrast, optical cluster studies, and in particular the MaxBCG 
group, frequently employ radii and masses scaled to the mean 
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Fig. 2. Scaled SZ signal measurements, Ysoo, binned by richness, A^200- The left-hand panel presents the results for the Johnston 
et al. (2007) M500 - N200 relation, the right-hand panel for the Rozo et al. (2009) relation. In each case, the red diamonds show 
the bin-average, redshift-scaled ?5oo calculated as the weighted mean of all individual measurements (e.g.. Fig. 1) in the bin, where 
the weights are taken from the estimated filter noise. The thick error bars show the corresponding uncertainty on the bin-average 
SZ signal, while the lighter error bars indicate the uncertainty found by bootstrap analysis; they are larger due to the presence of 
intrinsic scatter within the bins, most notable at high richness (see Fig. 4). The blue points represent the model prediction for each 
bin found by averaging, with the same weights as the data, the SZ signal expected from the Y^qo - M500 (Arnaud et al. 2010, STD 
case) and corresponding M500 - A^200 relations. The Planck measurements are little affected by choice of mass-richness relation, 
while the model points move significantly upward with the Rozo et al. (2009) mass calibration. Dashed lines in both panels show 
the best fit power-law to the Planck individual cluster data points (i.e., prior to binning, as shown in Fig. 1); the parameters for these 
fits are given in Table 2. 



matter density, rather than the critical density. For example, it 
is standard practice to refer to quantities measured within /?20()b, 
where the overdensity of 200 is defined with respect to the back- 
ground density (this corresponds to R(,o at z = and 7? 155 at 
Z = 1). For richness we will use the MaxBCG N200, defined 
as the number of red galaxies with L > 0.4L, within 7?2oob- 
Richness A^2()0 is the only quantity in this work defined relative 
to the mean background density. 

We characterize the SZ signal with the Compton-}; parameter 
integrated over a sphere of radius R^qq and expressed in arcmin^: 

5^500 - (crT/nisC^) ™ PdVID\{z), where Dp, denotes angular 
distance, cr-i is the Thomson cross-section, c the speed of light, 
me the electron rest mass and P - n^kT is the pressure, defined 
as the product of the electron number density and temperature, 
k being the Boltzmann constant. The use of this spherical, rather 
than cylindrical, quantity is possible because we adopt a tem- 
plate SZ profile when using the matched filter (discussed below). 
We bring our measurements to z = and a fiducial angular dis- 
tance assuming self-similar scaling in redshift. To this end, we 
introduce the intrinsic cluster quantity (an "absolute SZ signal 
strength") Y^qq = y5oo£'""^^(z)(I>A(z)/500Mpc)-, also expressed 
in arcmin^. 

2. Data Sets 

We base our study on Planck SZ measurements at the positions 
of clusters in the published MaxBCG cluster catalogue. 

2.1. The MaxBCG Optical Cluster Catalogue 

The MaxBCG catalogue (Koester et al. 2007b,a) is derived from 
Data Release 5 (DR5) of the Sloan Digital Sky Survey (York 
et al. 2000), covering an area of 7500 deg- in the Northern hemi- 



sphere. Galaxy cluster candidates were extracted by color, mag- 
nitude and a spatial filter centered on galaxies identified as the 
Brightest Cluster Galaxy (BCG). The catalogue provides posi- 
tion, redshift, richness and total luminosity for each candidate. 
In the following we will only use the richness A^200, defined as 
the number of red-sequence galaxies with L > 0.4L, and within a 
projected radius at which the cluster interior mean density equals 
200 times the mean background density at the redshift of the 
cluster (see Koester et al. (2007a) for details and the remark in 
Section 1.1). The catalogue consists of 13,823 galaxy clusters 
over the redshift range 0.1 < z < 0.3, with 90% purity and 85% 
completeness for 10 < A^200 < 190 as determined from simula- 
tions. 

A valuable characteristic for our study is the wide mass 
range spanned by the catalogue. Another is the fact that numer- 
ous authors have studied the catalogue, providing extensive in- 
formation on its properties. In particular, Sheldon et al. (2009) 
and Mandelbaum et al. (2008a) have published mass estimates 
from weak gravitational lensing analyses, which Johnston et al. 
(2007) and Rozo et al. (2009) use to construct mass-richness 
(M500 -N200) relations. We apply this relation, as outlined below, 
to adapt our SZ filter measurements for each individual cluster 
according to its given richness, A^200, as well as in our model 
predictions. 

In their discussion, Rozo et al. (2009) identify the dif- 
ferences between the Sheldon et al. (2009) and Mandelbaum 
et al. (2008a) mass estimates and the impact on the deduced 
mass-richness relation. They trace the systematically higher 
mass estimates of Mandelbaum et al. (2008a) to these authors' 
more detailed treatment of photometric redshift uncertainties 
(Mandelbaum et al. 2008b). Moreover, they note that Johnston 
et al. (2007), when employing the Sheldon et al. (2009) mea- 
surements, used an extended MaxBCG catalogue that includes 
objects with A^2oo < 10, where the catalogue is known to be in- 
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complete. These two effects lead Rozo et al. (2009) to propose a 
flatter mass-richness relation with higher normalisation than the 
original Johnston et al. (2007) result. In the following, we per- 
form our analysis with both relations; specifically, using the fit 
in Table 10 for the Mgoo - A^200 relation of Johnston et al. (2007), 
and Eqs. (4), (A20) and (A21) of Rozo et al. (2009). 

2.2. Planck Data 

We use the six HFI channel temperature maps (prior to CMB re- 
moval) provided by the DPC and whose characteristics are given 
in Planck HFI Core Team (2011b). These maps correspond to 
the observations of intensity in the first ten months of survey by 
Planck, still allowing complete sky coverage. Hence, they give 
us access to the entire SDSS survey area and complete MaxBCG 
catalogue. After masking bad pixels and contaminated regions 
(e.g., areas where an individual frequency map has a point source 
at > lOcr), we have Planck observations for 13,104 of the 13,823 
clusters in the MaxBCG catalogue. 

3. SZ Measurements 

We extract the SZ signal at the position of each MaxBCG clus- 
ter by applying a multi-frequency matched filter (Herranz et al. 
2002; Melin et al. 2006) to the six Planck temperature maps. 
The technique maximises the signal-to-noise of objects having 
the known frequency dependence of the thermal SZ effect and 
the expected angular profile. The filter returns the amplitude 
of the template, which we then convert into integrated SZ sig- 
nal, F500, within 7?5oo. It also returns an estimate of the local 
noise through the filter, ctq^^, due to instrumental noise and as- 
trophysical emissions. The same procedure is used in Planck 
Collaboration (201 If). We refer the reader to Melin et al. (2006, 
2010) for details. 

3.7. SZ Model Template 

For the filter's spatial template we adopt the empirical universal 
pressure profile of Arnaud et al. (2010), deduced from X-ray 
studies of the REXCESS cluster sample (Bohringer et al. 2007): 

P(r) oc (1) 

where the physical radius r is scaled to x - r/r^, with 
'"s - ^50o/c5oo- For the standard self-similar case [ST case in 
Appendix B of Arnaud et al. (2010)], C500 - 1.156 and the expo- 
nents are a = 1.0620,/? = 5.4807, y = 0.3292. The normaHsa- 
tion is arbitrary for purposes of the matched filter. The SZ signal 
being proportional to the gas pressure, we find the filter template 
by integrating along the line-of-sight and expressing the result 
in terms of projected angles: x - 919^. We truncate the filter at 
50500, containing more than 95% of the signal for the model. 

3.2. Application of the Filter 

We apply the matched filter to each cluster in the MaxBCG cat- 
alogue, using the mass-richness relation, M500 - A^200, to define 
R5Q0 and set the angular scale 6500 - R5qo/Da{z). The filter ef- 
fectively samples the cluster SZ signal along a cone out to a 
transverse angular radius of 5^?5oo, and returns the normalisation 
for the template. We apply a geometric factor based on the tem- 
plate SZ profile to convert the deduced total SZ signal along the 
cone to an equivalent Ksoo value, the SZ signal integrated within 



Table 1. Scaled Planck SZ signal measurements F500 binned by 
A^200 for the Rozo et al. (2009) mass-richness relation. Given Fgoo 
values are the measurement-noise weighted mean in the bin. The 
statistical uncertainty corresponds to the measurement-noise un- 
certainty on the weighted mean, while the total uncertainty ex- 
presses the standard deviation of the weighted mean from an 
ensemble of bootstrap samples. This table is plotted as the red 
diamonds and error bars in the right-hand panel of Figure 2. 



Moo 


?50o/(10"^arcmin^) 


Stat. Uncertainty 


Total Uncertainty 


10-13 


2.0 


±0.3 


±0.3 


14-17 


3.8 


±0.6 


±0.6 


18-24 


8.2 


±0.7 


±0.7 


25-32 


15 


±1 


±1 


33-43 


27 


±2 


±2 


44-58 


48 


±3 


±4 


59-77 


76 


±4 


±8 


78-104 


190 


±9 


±40 


>105 


300 


±20 


±80 



a sphere of physical radius R^oo- To account for the redshift range 
of the catalogue, we scale these measurements according to self- 
similar expectations to redshift z - and a fiducial angular 
distance of 500Mpc: F500 = F50o£"^^^(z)(£>A(z)/500Mpc)-. We 
accordingly adapt the estimated filter noise crg^^ to uncertainty 
CTgju^ on these scaled SZ signal measurements. The results of this 
procedure when using the Johnston et al. (2007) mass-richness 
relation are shown in Figure 1 . 

3.3. Systematic Effects 

As in the other four Planck SZ papers (Planck Collaboration 
201 ld,e,f,g), we have carried out various tests to ensure the ro- 
bustness of the Planck SZ measurements. They included inves- 
tigation of the cluster size-flux degeneracy, evaluation of the im- 
pact of the assumed pressure profile used for the Planck clus- 
ter detection, of beam-shape effects, color corrections, poten- 
tial contamination by point sources, as well as an overall er- 
ror budget estimation. We refer the reader to Sec. 6 of Planck 
Collaboration (201 Id) for an extensive description of this com- 
mon analysis. 

To complete this investigation in the present work, we re- 
peated our entire analysis, changing both the instrument beams 
and adopted SZ profiles. In the former instance, we varied the 
beams at all frequencies together to the extremes of their asso- 
ciated uncertainties as specified by the DPC (Planck HFI Core 
Team 201 lb). All beams were increased or all decreased in lock- 
step to maximize any effect. To investigate the profile, we re- 
extracted the SZ signal using a non-standard SZ signal-mass 
scaling, and separately for cool-core and morphologically dis- 
turbed SZ profiles (based on the work of Arnaud et al. (2010)). 
In all cases, the impact on the measurements was of order a few 
percent and thus negligibly impacts our conclusions. 

4. Results 

Our basic measurements are the set of individual scaled SZ sig- 
nal values Fsoo for each MaxBCG cluster, given as a function of 
richness A^200 in Figure 1 for the Johnston et al. (2007) mass cal- 
ibration. At high richness we can detect by eye a slight upturn 
of the points. Except for the most massive objects, however, the 
signal-to-noise of the individual measurements is small, in most 
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Fig. 3. Null test performed by randomising the angular positions of the clusters. The red diamonds show the bin-average, redshift- 
scaled measurements, F500, as reported in the left-hand side of Figure 2 with their corresponding measurement and bootstrap un- 
certainties; blue stars are the same model points. The green triangles present the bin-averages for the randomised catalogue with 
uncertainties given only by the SZ measurement errors. Results for the randomised catalogue are consistent with zero within their 
uncertainties. By comparison, the values for the real catalogue represent highly significant detections of the SZ signal in all richness 
bins. Left-hand panel: Results over the full richness range. Right-hand panel: Zoom into the region indicated by the dotted lines in 
the left-hand panel to highlight the low-richness end. 



cases well below unity. This is as expected given the masses of 
the clusters and the Planck noise levels. 

To extract the signal, we bin these F500 values by richness 
and calculate the bin averages as the noise-weighted mean of 
all individual / = 1, .., A^b measurements falling within the bin: 
<?500>b = {ZiY50o(i)/d-lji))/(Zil/d-lji)). We plot the re- 
sult as the red diamonds in Figure 2. The bold error bars repre- 
sent only the statistical uncertainty associated with the SZ signal 
measurements: crZ^ = y,- l/d-l (;) (in some cases the error bars 
are hidden by the size of the data point in the figure). The left- 
hand panel of the figure shows results using the Johnston et al. 
(2007) mass calibration, while the right-hand side gives results 
for the Rozo et al. (2009) mass calibration. The individual SZ 
signal measurements are not sensitive to this choice: the differ- 
ent calibrations do modify the adopted filter size, but the impact 
on the measured signal is small. 

We quantify the significance of the SZ detection using a null 
test: we perform an identical analysis on the MaxBCG cata- 
logue after first randomising the cluster angular positions within 
the SDSS DR5 footprint. In this analysis we are therefore at- 
tempting to measure SZ signal with the same set of filters, but 
now positioned randomly within the SDSS survey. The result is 
shown in Figure 3 by the green triangles, to be compared to the 
actual MaxBCG measurement given by the red diamonds. The 
left-hand panel presents the null test over the full richness range, 
while the right-hand panel affords an expanded view of the low 
mass end. The analysis on the randomised catalog remains con- 
sistent with zero (no detection) to within the SZ measurement 
uncertainty over the entire richness range. The actual measure- 
ments of the MaxBCG clusters, on the other hand, deviate by 
many cr from zero. We reject the null hypothesis in all bins at 
high significance. 

Figure 4 summarises our analysis of the uncertainty and in- 
trinsic scatter as a function of richness. In the left panel we show 
the uncertainty on the mean signal Fgoo in each bin, expressed 
as a fraction of Fsoq. The red solid red line traces the uncertainty 
on the mean signal due to just the measurement error, i.e., the 
noise level in the filter. The blue dashed line gives the uncer- 



tainty on the mean assuming that the measurements within a bin 
are Gaussian distributed about the mean with variance equal to 
the empirical in-bin variance. We show the relative uncertainty 
calculated from a bootstrap analysis of the entire catalogue as 
the dot-dashed, green curve. We perform our full analysis on 
10,000 bootstrap realisations from the actual catalogue and use 
the distribution of the resulting bin averages to find the relative 
uncertainty. The difference between the bootstrap and measure- 
ment uncertainties (red line) towards higher richness represents 
a detection of intrinsic scatter in those bins. At A^20o < 30, this 
difference is small and any intrinsic scatter is diflicult to distin- 
guish from the measurement errors. 

In the right-hand panel of Fig. 4 we show our estimate of 
the intrinsic scatter in the scaling relation as a function of rich- 
ness for A^200 > 30. This is expressed as a fraction of the mean, 
fsoo- The dot-dashed, blue line traces the empirical, or raw, dis- 
persion around the average signal of each bin. The three-dot- 
dashed, green line gives the dispersion corresponding to pure SZ 
measurement noise. To find the intrinsic scatter, we use the esti- 
mator: 

^ = ]^ Z (^500(0 - [?500]an.h) " " ^ <„(0 (2) 

where [?5oo]arith is the straight arithmetic mean in the bin. In the 
figure we plot 2b/(?5oo)b, with (?5oo)b being the weighted mean, 
as above. For this calculation we clip all outliers at > 5cr, where 
cr is the individual cluster SZ signal error. The final result, es- 
pecially at low richness, depends on the chosen clipping thresh- 
old. The scatter is not Gaussian, as the large fractional intrinsic 
scatter at low richness suggests. Below A^200 ~ 30, it becomes 
difficult to draw clear conclusions concerning the scatter, as can 
be appreciated by the fact that the bootstrap and pure SZ mea- 
surement uncertainties begin to overlap in the left-hand panel. 
For this reason, we only calculate the intrinsic scatter for the five 
highest richness bins in the right-hand panel. 

In conclusion, we detect a signal down to the lowest mass 
systems in the MaxBCG catalog with high statistical signifi- 
cance. This is the central result of our study. According to the 
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Fig. 4. Dispersion analysis. Left-hand panel: relative uncertainty on the mean versus richness. The relative uncertainty is expressed 
as a fraction of the bin-average redshift-scaled SZ signal: cr/Fjoo. The lower red curve corresponds to pure measurement uncertain- 
ties from the matched filter noise estimations; they are the solid error bars of Fig. 2. The upper blue curve traces the uncertainty on 
the mean assuming the points within a richness bin are normally distributed according to the observed in-bin dispersion. Bootstrap 
uncertainties are given as the middle green line, found as the dispersion in the mean Ysoq in each bin calculated over 10,000 boot- 
strap realisations of the entire MaxBCG catalogue. The numbers given in the legend indicate the number of objects in each richness 
bin. Right-hand panel: Fractional intrinsic scatter as a function of richness. The blue dot-dashed line (connecting the blue triangles) 
shows the raw dispersion in each richness bin, while the green dash-three-dotted line (connecting the green crosses) gives the cal- 
culated statistical dispersion from the measurement error on the scaled SZ signal Fsoq. The red dashed line with error bars is our 
estimation of the intrinsic scatter as a function of richness. For this calculation we have eliminated outliers in each bin at > 5cr, 
with o" — d-g^g„ for each cluster. We only calculate the intrinsic scatter at A^200 > 30, because at lower richness it becomes difficult to 
separate the intrinsic dispersion from the scatter due to pure measurement error. 



mass calibration from Johnston et al. (2007), we observe the SZ 
signal in objects of mass as low as M500 = (4 - 5) x lO'-' Mq. 



5. Discussion 

Figure 2 summarises the central results of our study. There 
are two notable aspects: firstly, we detect the SZ signal at 
high significance over the entire mass range; moreover, simple 
power laws adequately represent the observed scaling relations. 
Secondly, we see a discrepancy in the ?50{)-A^2(K) relation relative 
to expectations based on X-ray models and either the Johnston 
et al. (2007) or Rozo et al. (2009) mass calibrations. 
Fitting a power law of the form 



Da{z) 
500 Mpc/ 



(3) 



directly to the individual scaled measurements (e.g.. Fig. 1), we 
obtain the results summarised in Table 2. The Rozo et al. (2009) 
mass calibration assigns a larger mass to the clusters, increasing 
the filter scale and augmenting the measured SZ signal, which 
we see as the slightly higher normalisation. These fits are plotted 
as the dashed lines in Fig. 2. The power laws satisfactorily rep- 
resent the bin-average trends. The reduced = 1.16(13,104-2 
degrees-of-freedom) in both cases is poor; this reflects the pres- 
ence of the intrinsic scatter, also evident by the larger uncertain- 
ties on the fit from the bootstrap analysis. 

The blue stars in Figure 2 represent the predictions of a 
model based on the Kgoo-Msoo relation from Arnaud et al. (2010) 
and the Johnston et al. (2007) (left) or Rozo et al. (2009) (right) 
M500 - A^200 mean scaling relation. It assumes a self-similar 
I500 - M500 scaling relation (STD case) calibrated on X-ray 
observations of the REXCESS cluster sample (Bohringer et al. 
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Fig. 5. The Ysoo - N200 relation for the MCXC X-ray subsam- 
ple. Thick lines give the statistical errors, while the thin bars are 
the bootstrap uncertainties. We find that the MCXC X-ray sub- 
sample matches the model predictions much better than the full 
sample, which maintains a clear off'set relative to the model, as 
seen in in Fig. 2. 



2007). This calibration is also consistent with WMAP observa- 
tions (Melin et al. 2010) and with the Planck analysis (Planck 
Collaboration 2011f,g). In each bin we average the model pre- 
dictions in the same way as the Planck observations: we find the 
model bin-average redshift-scaled SZ signal as the inverse-error- 
weighted (pure SZ measurement error) average, assigning each 
cluster in the bin the same error as the actual observation of that 
object. Note that in the observation plane (?5oo, A^20o), the model 
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Table 2. The SZ signal-richness relation fit to a power law of the form F500 = Y20 (A^20o/20)" [see Eq. (3)] for the two mass-richness 
relations. The power law is fit directly to the individual SZ measurements (e.g.. Fig. 1). The columns labeled "Statistical" give the 
uncertainty on the parameters calculated from the SZ measurement errors alone, while those labeled "Bootstrap" give uncertainties 
found by fitting the power law to a set of bootstrap samples; the latter better represent the fuU uncertainty of the fits in the presence 
of intrinsic scatter. 



Mass-Richness Relation 


720/(10"' ai-cmin-) 


Statistical 


Bootstrap 


a 


Statistical 


Bootstrap 


Johnston et al. (2007) 


6.8 


±0.3 


±0.4 


2.07 


±0.03 


±0.07 


Rozo et al. (2009) 


7.4 


±0.3 


±0.4 


2.03 


±0.03 


±0.07 



(blue) points change with the mass calibration much more than 
the measurements. 

We see a clear discrepancy between the model and the 
Planck SZ measurements for both mass calibrations. In the case 
of the Johnston et al. (2007) mass calibration, the discrepancy 
manifests as a shift in normalisation that we can characterise by 
a 25% mass shift at given SZ signal: M — > 0.75M; the slope 
of the observed relation remains consistent with the self-similar 
prediction. The Rozo et al. (2009) mass calibration, on the other 
hand, flattens the mass-richness relation and predicts a shallower 
power law, as well as a higher normalisation; at A^200 = 50 there 
is a factor of 2 between the predicted and observed amplitudes. 

We now discuss some possible explanations for this discrep- 
ancy. Weak lensing mass estimates are difficult, and as we have 
seen there is an important difference in the two mass calibrations. 
Rozo et al. (2009), building on earlier work by Mandelbaum 
et al. (2008b), discuss some of the issues when measuring the 
weak-lensing signal for the MaxBCG catalogue. However, it 
seems unlikely that the weak-lensing mass calibration would be 
in error to the extent needed to explain the discrepancy seen in 
Figure 2. The discrepancy is in fact larger for the Rozo et al. 
(2009) result, which should be the more robust mass calibration. 

Our model predictions use a series of non-linear, mean rela- 
tions between observables which in reality have scatter that may 
also be non-Gaussian. The largest scatter is expected to be in 
the mass-richness relation. If the scatter is large enough, it could 
bias the predictions. We have investigated the effect of a 45% 
log-normal scatter in mass at fixed richness (e.g., Rozo et al. 
(2009)) and of a Poissonian distributio in richness at fixed mass. 
These are realistic expectations for the degree of scatter in the re- 
lations. The effect on the predicted, binned SZ signal is at most 
20%, not enough to explain the factor of two discrepancy we see. 

Contamination of the MaxBCG catalogue with a fraction, /, 
of objects that do not contribute an SZ signal (e.g., projection ef- 
fects in the optical) would bias the measured signal low by about 
1 - /. The level of contamination needed to explain the magni- 
tude of the discrepancy with the Rozo et al. (2009) calibration 
(/ 0.5) seems unlikely. The catalogue is estimated, instead, to 
be close to 90% pure for A^200 > 10. Moreover, contamination 
would also lower the weak-lensing mass calibration by about 
I -/, at given A^200- Since the predicted SZ signal scales as M^^^, 
the model SZ signal would drop by an even larger amount than 
the observed signal. 

To investigate this discrepancy further we analyse, in the 
same manner, a subsample of the MaxBCG clusters with X-ray 
data from the MCXC catalogue (Piffaretti et al. 2010). This rep- 
resents an X-ray detected subsample of the MaxBCG. The re- 
sults are given in Fig. 5 for the Rozo et al. (2009) mass calibra- 
tion and with our usual notation. We see that this X-ray subsam- 
ple, of 189 clusters, matches the model predictions much better. 
This argues that, at least for this subsample, the weak-lensing 
mass calibration is not significantly biased. The result hints at 



the presence of two populations: an X-ray under-luminous sam- 
ple, with a low SZ signal normalisation, seen only in the opti- 
cally selected catalogue, and an X-ray normal sample matching 
X-ray model predictions. We emphasize, however, that the se- 
lection function of this X-ray subsample is not known, so this 
interpretation is provisional. 

Splitting the catalogue according to the luminosity of the 
BCG lends support to the presence of two populations. In each 
bin, we divide the catalogue into a BCG-dominated sample, 
where the fraction of the cluster luminosity contributed by the 
BCG is larger than the average for that bin, and its complement 
sample. The BCG-dominated sample has a notably higher nor- 
malisation, closer to the predicted relation, than the complement 
sample. 

We also compare our results to the X-ray results from Rykoff 
et al. (2008a) who stacked ROSAT photons around MaxBCG 
clusters according to richness. As with our SZ observations, their 
individual X-ray fluxes had low signal-to-noise, but they ex- 
tracted mean luminosities from each image stack (Rykoff et al. 
2008b). They report luminosities, L200, over the 0.1 - 2.4keV 
band and within R200 for each A^200 richness bin. Their analysis 
revealed a discrepancy between the observed mean luminosities 
and the X-ray model predictions, using the Johnston et al. (2007) 
mass calibration. 

To compare we re-binned into the same richness bins as 
Rykoff et al. (2008b), calculating the new, bin-average, redshift- 
scaled Fsoq. We also convert their luminosities to L500 using the 
X-ray profile adopted in Arnaud et al. (2010); the conversion 
factor is 0.91. In addition, we apply the self-similar redshift lu- 
minosity scaling of E^''^^{z = 0.25) to bring the Rykoff et al. 
(2008a) measurements to equivalent z - values from the val- 
ues at their median redshift, z - 0.25. The resulting points are 
shown in Fig. 6. 

The model line in the figure is calculated from the z = 0, X- 
ray luminosities using the X-ray based scaling laws in Amaud 
et al. (2010). In this plane the model matches the observations 
well, demonstrating consistency between the SZ and X-ray ob- 
servations. Remarkably, the ICM quantities remain in agreement 
with the model despite the individual discrepancies (SZ and X- 
ray luminosity) with richness. 

The intrinsic scatter in the scaling relation, given in Figure 
4, starts at about 60% and rises to over 100% at A^200 ~ 30. 
This was calculated by clipping all outliers at > 5cr; the re- 
sult depends on the choice of clipping threshold, indicative of a 
non-Gaussian distribution. This dispersion should be compared 
to the estimated log-normal scatter in the mass-richness relation 
of (45!ojg)% fo™'^ by ^1- (2009). Assuming that the 

dispersion in the SZ signal-mass relation is much smaller, we 
would expect a dispersion of order 75%, not far from what we 
find and within the uncertainties. Such large fractional disper- 
sion implies a non-Gaussian distribution skewed toward high SZ 
signal values, particularly at low richness. 
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6. Conclusions 

We have measured with high significance the mean SZ signal for 
MaxBCG clusters binned by richness, even the poorest systems. 
The observed SZ signal-richness relation, based on 13,104 of 
the MaxBCG clusters observed by Planck, is well represented 
by a power law. This adds another scaling relation to the list 
of such relations known to exist among cluster properties and 
that present important constraints on cluster and galaxy evolu- 
tion models. 

The observed relation has a significantly lower amplitude 
than predicted by X-ray models coupled with the mass-richness 
relation from weak-lensing observations. The origin of this dis- 
crepancy remains unclear. Bias in the weak-lensing mass mea- 
surements and/or a high contamination of the catalogue are po- 
tential explanations; another is the presence of a wide range of 
ICM properties at fixed richness, of which only the more X-ray 
luminous are readily found in X-ray samples used to establish 
the X-ray model. The better agreement of the model with a sub- 
sample of the MaxBCG catalogue with X-ray observations sup- 
ports this conjecture, although the unknown selection function 
of the X-ray subsample renders it provisional. Remarkably, the 
relation between mean SZ signal and mean X-ray luminosity for 
the entire catalogue does conform to model predictions despite 
discrepant SZ signal-richness and X-ray luminosity-richness re- 
lations; properties of the gas halo appear more stably related than 
either to richness. 

We find large intrinsic scatter in the SZ signal-richness re- 
lation, although consistent with the major contribution aris- 
ing from scatter in the mass-richness relation. The uncertain- 
ties, however, are important. Such large scatter implies a non- 
Gaussian distribution of SZ signal at given richness, skewed to- 
wards higher signal strengths. This is consistent with the idea 
of a wide range of ICM properties at fixed richness, with X-ray 
detected objects preferentially at the high SZ signal end. 

If the lower normalisation of the Fsoo - A^200 relation, and by 
consequence the ?5oo - A^soo relation, is confirmed as a physical 
and representative property of the cluster population, as opposed 
to selection or other observational effects, then this would have 
important implications for our understanding of clusters. Among 
others, it would imply a lower amplitude for both the diffuse SZ 
power spectrum and for number counts of SZ-detected clusters. 
Predictions for these two quantities depend on the ^500 - M500 
relation. The amplitude of the SZ power spectrum varies as the 
square of the normalisation, while the counts depend on it expo- 
nentially. In both instances, this relation represents a significant 
theoretical uncertainty plaguing models. 

Our study of the SZ signal-richness relation is a step to- 
wards reducing this uncertainty, and it presents a new cluster 
scaling relation as a useful constraint for theories of cluster and 
galaxy evolution. Concerning the latter, we find no obvious sign 
of an abrupt change in the ICM properties of optically selected 
clusters over a wide range of richness, hence mass, as might 
be expected from strong feedback models. Future research with 
Planck will extend this work to other catalogues and a greater 
redshift range. 
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Fig. 6. Comparison of our bin-average redshift-scaled SZ sig- 
nal measurements with the mean X-ray luminosities found by 
Rykoff et al. (2008a). For this comparison, we re-bin into the 
same bins as Rykoff et al. (2008a) and plot the results as the red 
diamonds with error bars. The X-ray luminosities are brought 
to equivalent z - values using the self-similar scaling of 
E^^^^iz - 0.25) applied at the quoted z = 0.25 median red- 
shift. The dashed blue line shows the predictions of the X-ray 
model. Our notation for the eiTor bars follows previous figures. 
The numbers in the figure indicate the number of clusters in each 
bin. 
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